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Set Books 


D. L. Kreider, R. G. Kuller, D. R. Ostberg and F. W. Perkins, An Intro- 
duction to Linear Analysis (Addison-Wesley, 1966). 
E. D. Nering, Linear Algebra and Matrix Theory (John Wiley, 1970). 


It is essential to have these books; the course is based on them and will 
not make sense without them. 


Conventions 


Before working through this correspondence text make sure you have 
read An Introduction to the Bridging Material and A Guide to the Linear 
Mathematics Course. 


The set books are referred to as: 


K for An Introduction to Linear Analysis 
N for Linear Algebra and Matrix Theory 


Note 


This bridging material is not based on the set books. It has been written 
especially for the benefit of students who have taken the Mathematics 
Foundation Course M101 (The Open University Press, 1978) 

References to this foundation course take the form M101 Block V Unit 2 
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21 TAYLOR POLYNOMIALS AND 
APPROXIMATIONS TO FUNCTIONS 


21.1 Introduction 


You will remember from M101 Block III Unit 4 p18 that the n* Taylor 
polynomial for a function f, provided f can be differentiated n times, is 


1 
FO) ESO x + FSO)? e fo 
You will also have seen in TV14 how to use the Taylor polynomials for 


the sine function to calculate approximate values of sin x: computer 
graphics demonstrated that the fifth Taylor polynomial for sine 


does not differ very much from sine x over most of the interval 


- 5, 5 . Because of the application of Taylor polynomials calculating 
approximate values of functions, the n'^ Taylor polynomial of a function 
f is called the »'^ Taylor approximation to f in M201 and elsewhere. 
(You may find it worth while to watch M101 TV14 again as you will 
need this work before reading M201 Units 19 and 21). 


The important information that we must have before we can use Taylor 
polynomials to calculate (approximate) values of sine, or of other 
functions, and have any confidence in our answer, is to know how 
accurately we can expect a Taylor polynomial to approximate the 
function. The difficulty is precisely the same as the one you met in M101 
Block II Unit 3 in discussing error bounds; in practice, if we are using 
the Taylor polynomial to calculate an approximate value for a function 
such as sin x; all we know is the approximate value, not the true value, 
so we don’t know the size of the error! The way out of this difficulty is 
given by Taylor’s theorem, which tells us how to calculate our upper 
limit, or upper bound for the error. This upper bound will usually over- 
estimate the actual error, but it will never the less be very useful. 


2.4.3. Taylor's Theorem (Simple Case) 


The first Taylor approximation (first Taylor polynomial) of a function f 
is f(0)--xf'(0) Compare this with the second approximation 
S (0) + xf'(0) + 3x?f"(0). The difference is 1x?/"(0) and this suggests that 
an approximation to the error in the value given by the first formula is 
3x?f"(0). Obviously this can't be the upper bound for the error that we 
are seeking, but in fact Taylor’s theorem tells us that the upper bound 
looks rather like it. The precise statement is as follows. 


Suppose f can be differentiated twice. Then 
J(=) = £(0) + xf'0) + C.) 

where the correction term C,(x) satisfies 
[Ci (x)| x 48x? 


provided |f"(r)| < B for all te [0, x] (or te [x, 0] if x <0) so if we 
replace f"(0) in the quadratic term 3x?f"(0) by a number B which is 
bigger than the value of f"(t) over the whole interval [0, x], we get a 
bound for the error in the first Taylor approximation. 
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We shall not prove this result here: you will find a proof of the general 
case on page K667 if you are interested, but a knowledge of the proof is 
not required for M201. 


Example 1 


To illustrate Taylor's theorem, let us apply it to the sine function. In this 
case, the first Taylor approximation is 


sin x œ sin 0 + x sin’ (0) 
=0+x 
and Taylor's theorem tells us that the error — C, (x) satisfies 
|C:69)] < $8x? 
provided |sin" (t) < B for all t e [0, x]. 


Since we know sin" (t) 2 —sin t and |sin t| < 1 for all t we can take 
B — 1. Thus Taylor's theorem tells us that 


[sin x — x| < 4x? 
For example, if x = 45, we can deduce sin 75 = 0.100 with error at most 
2 
35)" = 0.005 


Notice that the true value of sin T to five decimal places is 0.09983 so 


that the actual error is only about one twentieth of that guaranteed by 
Taylor's theorem. However, we do have for certain a maximum value for 
the error, and we got it for very little work. Indeed, as we shall see 
shortly, this error estimate could have been improved considerably with 
only a little more care. 


Example 2 
big ~ 


Show that the error in taking tan = é 0.52 (the value of the first 


Taylor approximation to tan x) is less than 


nf = 0.21 
If J (x) = tan x, 
then — f'(x) = sec? x 

I(x) = 2 sec? x tan x 


so the first Taylor approximation to tan x is 0 + xf'(0) = x. 


2 
The error for x= is at most T provided |f"()| < B for 
T 
t F2 
€ le 6 


Since sec t and tan t both increase as t increases on |0, 


d the largest 


value that f"(t) can take on the interval b. d is 


2V 1 8 
2 sec? * tan 7 2x ( ) x = à 
6g JA aA 
1 2 x? 
Hence the error is at most 2: —2--: 7; = ~ 021. 


8 m 
23/360 21/3 
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2.13 Taylor's Theorem (General Case) 


The simple case of Taylor's theorem can be generalized in two ways. 
First, we can obtain an upper bound for the error in using the n" 
Taylor approximation. 


Theorem 
Suppose f is differentiable (n + 1) times. Then 


FO) = £00) + x0) 5 $0) + c.) 


where the correction term C,(x) satisfies 


xt 
(n+ 1)! p 
provided | f@*(t)| < B for all re [0, x] (or t e [x, 0] if x « 0). 


IC x) < 


Example 3 


Use Taylor's theorem to find an upper bound for the error in using the 
second Taylor approximation to calculate sin 0.1. 


Solution 


Since sin’ (0) = cos 0 = 1, 
sin" (0) = —sin 0 = 0 


The second Taylor approximation to sine is the same as the first namely 
Sin x = x. 


We can write sin x = x + C;(x) where by Taylor's theorem C ;(x) satisfies 


ICa(x)| <| 57 B 


provided |sin” (t)| < B for all t e [0, x]. 


Now sin" (t) = —cos t and [cos t| < 1 for all t, so 


[C369] < 


x? 
3! 


In particular, sin 0.1 = 0.1 + C4(0.1) 


3 
where |C;(0.1)| < e». 0.00017 


This answer shows a considerable improvement on the error bound 
obtained in Example 1, although the approximate value of sin (0.1) is the 
same in each case. 


The other way in which we can generalize Taylor's theorem is so that 
we can use it not only when we can easily work out f and its derivatives 
at x = 0, but also where x = o is the point where we know something 
about the function. This is just like the situation in M101 Block III 
Section 4.4 where we considered Taylor series about a general point. 
The result that we obtain is 


fe) = f) + Sla) — o) A fl) — a)? 
1 


ed f" =a) +e " 


Lax — a)” Ts 


2.4.4 "Taylor's Theorem 
If f is differentiable (n + 1) times then 
ft) = f) (— 9/6): FEM" pomay + ca) 


where C,(x) satisfies 


yy 
(x-aytt a 


le) < 


provided | f*"(t)| < B for all t e [a, x] (or [x, a] if x <a). 


Example 4 

Use the second Taylor approximation at x = 1 to calculate log, 1.1. 
Find an upper bound for the error, using Taylor's theorem. 

Solution 

If f(x) = log, x 

then. f'(x)2x^! 


f'e)- =x"? 
fi) 2 


The second Taylor approximation at x — 1 is 
Se) = SU) + (x — DFC) + Hoe — 1700) 
-0-4(x—1)-3x- 1)? 
so f (11) = 0.1 — 40.1)? = 0.9950. 
By Taylor's theorem, the correction term C (x) satisfies 
[C6] < [a(x — 1)°B] 
provided | f”(t)| < B for all te [1, x]. 
Here, x = 1.1 and f"(t) = 217? <2 for te [1, 1.1]. 


Hence the error in the estimate 1og,1.1— 0.095 is less than 
$(0.1)*.2 < 0.0004. 


You can now read M201 Unit 19, Section 3.4 on Taylor Approximation. 


22 PARTIAL DIFFERENTIATION 


2..1 Introduction 


In the case of a function of a single variable, we are able to measure 
the rate of change of the value of the function relative to a change in 
the value of the variable by finding the derivative of the function. In 
the case of a function of two variables, we can seek similarly to measure 
the rate of change of the value of the function relative to a change in the 
value of either of the variables. The new concept we need is that of a 
partial derivative. 


It is easy to get an intuitive grasp of the notation of a partial derivative. 
Imagine yourself standing on a hillside at a point where two roads cross, 
one road running east-west, the other north-south. One road may perhaps 
have a steep gradient, the other less so. The surface of the hill represents 
the graph ofa function of two variables, Northings and Eastings. (Height 
depends upon both the latitude and the longitude of the point in question). 
Speaking roughly, the two partial derivatives of the function ‘height’ at 
the crossroads are measured by the slopes of the two roads. If the 
crossroads happened to be at the top of the hill then each of the slopes 
would be zero at that point, but elsewhere they will differ. 


We want to make this intuitive idca more precise. The geometric 
example which follows shows the way we must go and helps you to 
understand the subsequent definitions. But, if you find it hard to visualize 
three-dimensional figures, do not spend a lot of time trying to understand 
the example. 


Example 


Consider the surface representing the function 


z-F(sy)-1-(Q y) 


(x y)E E x R, x? + y? < 1) 


We always take the positive square-root. Then, if we consider z to be 
‘height’ and x and y to be ‘Eastings’ and ‘Northings’, the graph of 
the function F corresponds to the hillside of our intuitive discussion. 


Since z = ,/1 — (x?  y?), it follows that x^ + y? + z? = 1 so that all 
points of the graph lie on the surface of a unit sphere with centre at the 
origin. 

Choosing z to be always positive means that its graph is in fact a 
hemisphere, looking a little like a hill. 

On this hill we choose, for illustration, the point Q with co-ordinates 


(3. 0, Y. Jana we take two roads through Q given by the intersection 
of the hemisphere with the vertical planes through Q parallel to the x 
and y axes. 


The whole of the east-west road lies in the plane y=0. 


If we were to cut the hemisphere with the plane y = 0 through Q. and 
then look along the y-axis, we would see the semi-circle shown in red 
in the following diagram: 


The plane y - 0 


The plane y = 0 


The semi-circle is the intersection of the plane y = 0 with the hemisphere 
z=,/1— (X? + y?) and so its equation, in the plane y = 0, is 


z= /1—x?, xe[-1, +1]. 


Defining a function f, by 


Filx) = /1 = x3, (xe[-L +1] 


10 
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we can calculate the slope of the road at Q from the derivative 


file) = Fs. 


1 ‘ 
When x=5, this takes the value M which is the slope of the 


V 


The north-south road similarly lies at the intersection of the plane x 2 4 


with the hemisphere z = ./1 — (x? + y?). 


We thus have 

z= SIEF V (ve 
Defining a function f; by 

fly) = /i- y. (ve -£ ex 


we calculate the slope from the derivative 


easterly road through Q. 


Aly)= g 


When y — 0 this takes the value 0, which is the slope of the northerly 
road through Q. : 


22.2 Definition of Partial Derivatives 


The intuitive idea of the previous example points the way to a more 
general definition. 


Consider a point P = (a, b, c) on the hemisphere, so that c > 0 and 
a? + b? + c? = 1. Through P we have the two planes y = 6 and x =a 
and we may seek the slope at P of the semi-circles in which each plane 
meets the hemisphere. 


Keeping y fixed at the constant value b, we obtain a function of just 
one variable: 


Silx) = 1- (x? + b?), xe[-/1 —b*, - /1— b7] 


This function has a graph which is a semi-circle in the plane y — b and 
the slope is given by'the derivative 


Similarly, if we keep x fixed at the constant value a, we get another 
function of a single variable: 


fix) = 


fax) = V1- (a + y?) ye[- /1- à, tI-— a] 


The graph of this function is a semicircle in the plane x — a and its 
slope is given by the derivative 


Domain of f, 


Domain of F 


By defining f;'(x) for each value of be [—1, +1], we may construct a 
new function of two variables: 


Zx 24.2 
JICQI («)e():x +y? «1 


Likewise, defining a function f;'(y) for each ae [—1, +1], we construct 
another function of two variables: 


Fi'(x, y)= 


F(x, y)= VEET (x, ye {(x, y): x? + y? < 1}. 


F,' and F, are the two partial derived functions of F. 


Exercise 


Calculate F,'(5, 0) and F,‘(3, 0) and compare your answers with the 
slopes obtained in the example of the previous section. 


Solution ; i 
F/8,0) 2 ———3 os 
JÀi-(y-o) Ji-i V3 

"(1 = -0 209. 

"= aro an” 


These values agree with those previously obtained for the slopes 


at Q. 


Remembering that all we do is to keep each variable in turn fixed whilst 
we differentiate with respect to the other variable, we can define the 
partial derivatives of a general function 


F:Rx ROR. 
The partial derivative of F with respect to the first variable, x, at the point 
(x, y) is 
F(x + h y) — F(x, y) 


F,'(x, y) = lim EMIL RE 
h-0 1 


The partial derivative of F with respect to the second variable, y, at the 
point (x, y) is 


3 . F(x, y +k) - F(x, 
Fs y) = tim FY + 16839) 


Example 1 

Consider the function 
G: (x, y) —2xy + x?, (x, y)E R x R. 

Treating y as a constant, we may differentiate with respect to x and get 
Gy'(x, y) 2 2y + 2x. 


Similarly, treating x as a constant, we may differentiate with respect to y 
and get 


Ga'(x, y) = 2x. 


Exercise 


Verify the results of the above example by working directly from the 
limit definitions. 


Solution 
, 2y(x + h) + (x + h}? — Qxy + x? 
aiig (EEN EEN tees =) 
hoo 
(= + 2xh+ “) 
Sim (=e 
ho h 
= 2y + 2x 
2 2 
Sess (2e +k)+ x (2xy +x )) 
k=0 
«e 
k-0 k 
- 2x. 
Exercise 


Each of the functions defined below has domain R x R. Find all their 
partial derivatives. 

() Fe y)=x?+y’. 

(ii) G(x, y) = x exp(xy). 

(iii) H(x, y) =x sin(x + y). 

(i) P(x, y) = x4 + y* — 4x?y). 


Solution 
Q FU) =2x, 

Fy, y) = 2y. 

(il) Gy'Gx, y)  exp(xy) + xy exp(xy), 
G2'(x, y) = x? exp(xy) 

(iii) H,'(x, y) = sin(x + y) + x cos(x + y), 
H'(x, y) = x cos{x + y). 

(iv) Pi'(x y) = x — 8xy%, 
P,'(x, y) = 4? — 12x?y?. 


Notation 


You will meet various alternative notations for partial derivatives. 
Thus F, is often used where we have F,’, and F, for Fz’. 


The commonest notation of all is 


oF 
x for F,'(x, y) and 


OF 
y for F2'(x, y). 


d ; -— 
This is reminiscent of the use ord for the ordinary derivative /'(r). 


However, be very wary of jumping to conclusions: it is not generally true 


that 2r i is the same as 7: whatever the notation may suggest. 
xX 


